Heterovalent interlayers and interface states: an ab initio study of 
GaAs/Si/GaAs (110) and (100) heterostructures 
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We have investigated ab initio the existence of localized states and resonances in abrupt 
GaAs/Si/GaAs (110)- and (lOO)-oriented heterostructures incorporating 1 or 2 monolayers (MLs) 
of Si, as well as in the fully developed Si/GaAs (110) heterojunction. In (lOO)-oriented structures, 
we find both valence- and conduction-band related near-band edge states localized at the Si/GaAs 
interface. In the (110) systems, instead, interface states occur deeper in the valence band; the high- 
est valence-related resonances being about 1 eV below the GaAs valence-band maximum. Using 
their characteristic bonding properties and atomic character, we are able to follow the evolution of 
the localized states and resonances from the fully developed Si/GaAs binary junction to the ternary 
GaAs/Si/GaAs (110) systems incorporating 2 or 1 ML of Si. This approach also allows us to show 
the link between the interface states of the (110) and (100) systems. Finally, the conditions for the 
existence of localized states at the Si/GaAs (110) interface are discussed based on a Koster-Slater 
model developed for the interface-state problem. 

PACS numbers: 73.20.-r, 73.20.At, 73.40.Ty 



I. INTRODUCTION 

In recent years, a large amount of theoretical and ex- 
perimental work has been concentrated on heterostruc- 
tures composed of ultrathin layers of a given group-IV 
material deposited at III-V/III-V or II-VI/II-VI homo- 
or hetero- junctions^ Interest in these systems has been 
stimulated by the possibility of tuning the band disconti- 
nuity in heterojunctions, or inducing a band discontinu- 
ity in homo junctions, thus allowing more flexibility in de- 
signing novel devices. 2 In particular, model calculations, 3 
supported by ab initio computations, 4 have indicated 
that Si inserted as a bilayer at a GaAs/AlAs (100) het- 
erojunction or at a GaAs/GaAs (100) homojunction, can 
act as a microscopic capacitor and induce substantial 
changes in the band discontinuities. Although the in- 
tcrlayer morphology is still difficult to controljMi sim- 
ilar band-offset modifications have also been reported 
experimentallyi&Silfi An issue, however, which has re- 
mained largely unexplored concerning such engineered 
interfaces is that of the electronic states induced by the 
interlayer. For any practical application, it is indeed im- 
portant to know whether the interlayer will induce lo- 
calized interface states within the bandgap that may de- 
grade the transport and optical properties of the junc- 
tion. 

Theoretical investigation of electronic states local- 
ized at interfaces between heterovalent materials has 
a long history beginning with the pioneering work of 
Baraff et al. 11 on the Ge/GaAs (100) interface. Since 



then, several different heterojunctions have been in- 
vestigated with special emphasis on the interface-state 
problemiiiii 2 ^ 3 ^ 4 * 1 ^* 1 ^ 1 ^ However, to our knowledge, no 
theoretical study of the interface states has been per- 
formed so far on the Si/GaAs heterojunction. Moreover, 
a complete understanding of the mechanisms that con- 
trol the formation of interface states and the conditions 
of their existence is still lacking. As far as the prob- 
lem of electronic states induced by heterovalent interlay- 
ers is concerned, to our knowledge, studies in this area 
have been limited to semiempirical tight-binding investi- 
gations on Ge interlayers in GaAs, and of Si and GaAs 
interlayers in II- VI semiconductorsi 1 ^* 1 ^* 1 ^ Tight-binding 
calculations, however, neglect charge effects associated 
with the heterovalent nature of the interlayerji 4 * 1 ^ and 
have yielded conflicting results (depending on the pa- 
rameters used) on the presence of localized states within 
the gap of polar systems, such as the Ge/GaAs (100) 
heterojunction. 13,14 

In the present work, we systematically investigate ab 
initio the localized states and resonances induced by thin 
layers of Si in GaAs (110) and (100) homojunctions, as 
well as in the fully developed Si/GaAs (110) heterojunc- 
tion. Our results reveal interesting connections between 
the formation of localized states and resonances in het- 
erostructures possessing different orientations and for in- 
terlayers of different thicknesses. A model is also pre- 
sented, which provides insight into the mechanisms of 
interface-state formation; it explicitly links the existence 
of the heterojunction interface states to some essential 
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bulk band-structure features of the constituent materials 
and some interface-specific bonding parameters^ Since 
the abrupt polar Si/GaAs (100) interface is charged, 
and therefore gives rise to a macroscopic electric field 
which complicates the study — and renders the fully de- 
veloped, abrupt (100) heterojunction thermodynamically 
unstable^ — we start by examining the interface-state 
spectrum of the non-polar Si/GaAs (110) heterojunction. 
We then study the electronic states induced by 1 and 
2 MLs of Si in the GaAs (110) homojunction, and by 
2 Si MLs in the GaAs (100) homojunction. We discuss 
the relationship between these states and their interface 
counterparts, and the conditions that govern the exis- 
tence of localized interface states. 



II. AB INITIO METHOD 

The interface states and interlayer-induced states are 
studied assuming coherently strained Si layers on a GaAs 
substrate. The theoretical framework is density func- 
tional theory (DFT) in the local density approximation 
(LDA), using the Ceperley- Alder exchange-correlation 
functional^ We use scalar relativistic Troullier-Martins 
pseudopotentials, 20 in the Kleinman-Bylander non-local 
form. 21 The valence electronic states are expanded on 
a plane-wave basis set, using the supercell technique 
to model the heterojunctions and heterstructures. We 
employ supercells including 12 monolayers (24 atoms) 
for the (llO)-oriented structures and 16 monolayers (16 
atoms) for the (100) structures. The supercell calcula- 
tions are performed with a kinetic-energy cutoff of 20 Ry 
and a (6,6,2) Monkhorst-Pack fc-point meshiSS For the 
bulk calculations we use a 40-Ry cutoff and a (6, 6, 6) 
mesh. The calculated equilibrium lattice constants of 
bulk GaAs and Si are 10.47 a.u. and 10.24 a.u., re- 
spectively. The resulting lattice mismatch is ~ 2%, to 
be compared with the experimental value of ~ 3. 5%. 23 
The theoretical values of the elastic constants of Si are 
given m Table |U they compare well with the experimen- 
tal values (also reported in the table) , and with previous 
calculations^ Using the calculated values of the elas- 
tic constants, and based on macroscopic elasticity the- 
ory (MET)^ the theoretical lattice constant of Si on 
GaAs along the growth direction is 10.13 a.u. [10.08 
a.u.] for the (llO)-oriented [(lOO)-oriented] structures. 
In the present work, for the starting (unrelaxed) Si/GaAs 
heterojunction and GaAs/Si/GaAs heterostructures, we 
have used the MET result for the Si-Si interlayer spacing, 
and, for the Si-GaAs interplanar distance at the inter- 
face, we have used the average between the MET inter- 
layer spacings in the two bulks. 

The supercell electronic states have been classified ac- 
cording to their degree of localization at the interface 
(or at the interlayer location) and their atomic-like na- 
ture. We have considered for the interface-state spectrum 
only those states whose probability density has a much 
higher value at the inversion layer than in the bulk re- 



TABLE I: Theoretical equilibrium lattice parameter and elas- 
tic constants of Si used in this work. Experimental values are 
reported in square brackets (from Ref. l23l) . 

ao (a.u.) en (Mbar) C12 (Mbar) C44 (Mbar) 
Si 10.24 [10.27] 1.62 [1.67] 0.56 [0.65] 0.73 [0.80] 



gions. Such states may be either truly localized states, 
which decay exponentially in the bulk materials, or reso- 
nances. Localized states are normally found in common 
gaps of the projected band structures (PBS's) of the bulk 
materials^ while resonances are degenerate with Bloch 
states of one or both bulks. It is therefore convenient to 
represent the interface states on a common energy scale 
with the PBS's of the bulk materials. The alignment of 
the bulk band structures was performed using the calcu- 
lated valence-band offset (VBO) obtained with the tech- 
nique outlined in Ref. |23- The supercell interface-state 
spectrum has then been aligned with respect to the top 
valence-band edge of one of the bulk materials using the 
calculated shift of the macroscopic-averag o 27 ' 28 electro- 
static potential in the corresponding bulk region of the 
supercell relative to the mean value of the electrostatic 
potential in the supercell. 

In the case of the (llO)-oriented structures, atomic re- 
laxation was found to have a relatively minor influence 
on the interface-state spectra. The choice of the average 
bulk Si and GaAs MET spacings for the Si-GaAs inter- 
layer distance is a good starting point for these struc- 
tures: the error with respect to the relaxed interlayer 
distance is less than 1%. Atomic relaxation has nor- 
mally two effects. It changes: (i) the band offset, i.e., 
the relative position of the bulk PBS's and (ii) the rel- 
ative position of the interface states with respect to the 
bulk PBS. The first effect vanishes (by symmetry) for the 
GaAs/Si/GaAs (110) heterostructures, and is negligible 
for the Si/GaAs (110) heterojunction (AV rs 0.02 eV). 
The deviations (ii), instead, can be as large as 0.1 eV for 
the Si/GaAs (110) heterojunction. However, this is of 
the same order of magnitude as the many-body and spin- 
orbit effects^ not included in our calculations. There- 
fore, for the (110) structures we will consider in the fol- 
lowing the ideal unrelaxed geometries only. 

Larger relaxation effects were found, instead, in 
the (lOO)-oriented GaAs/2 SiMLs/GaAs heterostruc- 
ture. The As-Si (Ga-Si) interplanar distance increases 
(decreases) by 7.5% (0.5%) and the Si-Si interlayer spac- 
ing decreases by 2.4%, with respect to the unrelaxed 
MET geometry. This results in changes of several tenths 
of an eV in the energy position of the localized states 
and resonances of the GaAs/2 SiMLs/GaAs (100) het- 
erojunction. In the following, all results presented for 
the GaAs/2 SiMLs/GaAs (100) system will refer to the 
fully relaxed geometry. 

We have also checked the effect of the supercell size 
using a 32-atom cell for the Si/GaAs (110) heterojunc- 
tion. We found variations in the energy position of the 
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electronic states of at most 0.05 eV. We can therefore 
conclude that the energy position of the electronic states 
in our (100) and (110) heterostructures is given with an 
overall uncertainty of the order of 0.1 eV. 



III. SI/GAAS (110) INTERFACE STATES 

Following the procedure outlined in the previous sec- 
tion, we have calculated the interface states of the 
Si/GaAs (110) heterojunction. The result is displayed in 
Fig. ^ where we have represented the dispersion of the 
states along the principal lines of the irreducible wedge 
of the two-dimensional Brillouin zone (2DBZ), together 
with the two bulk PBS's. The calculated VBO of the 
Si/GaAs (110) heterojunction (without spin-orbit and 
self-energy corrections) used for the alignment of the bulk 
PBS's is 0.14 eV (Si higher). The high-symmetry points 
of the_2DBZ are T = (0,0), X = (V%0), M = (\/2,l), 
and X' = (0,1), in units of ir/a (a is the GaAs lattice 
constant). 

The isolated (110) interface has a a yz reflection 
symmetry^ (z is the growth direction and y is the di- 
rection perpendicular to the Si-Si or Ga-As nearest- 
neighbor zigzag chains in the (110) atomic planes). How- 
ever, since we have two interfaces per supercell an ad- 
ditional (fictitious) reflection symmetry induces non- 
degenerate pairs of symmetric and antisymmetric (by re- 
flection) states in the quasi 2DBZ of the supercell (the 
symmetry group of the supercell is C^u). With the su- 
percell employed, though, the energy difference between 
the above pairs is about 50 meV throughout the BZ, thus 
indicating a small interaction between the two interfaces. 

Beginning from the lowest energies, we find a complete 
band of localized states below (or at the edge of) the bulk 
PBS of the s valence band of GaAs; the probability den- 
sity of the state at the X ' point of the 2DBZ is shown 
in Fig.[2Ia). These states (Sas in Fig. ^) derive from the 
s states of the As atoms adjacent to the interface. The 
localization of these states results from the more attrac- 
tive ionic potential of Si relative to Ga, which lowers the 
energy of the As-s states at the interface with respect 
to those in bulk GaAs. This is consistent with the rel- 
ative strengths of the atomic potentials, which obey the 
inequality Ga < Si < As, and leads to the following ex- 
pected attractiveness of the potential in the interfacial 
bonding regions: Si-Ga < Si-Si < Ga-As < Si-As. We 
note that similar localized Sa s states were reported by 
Pickett et alM for the Ge/GaAs (110) interface. 

Moving upwards in energy we find a second band of 
states which have a maximum probability density on the 
first Si bilayer from the interface (Ssi states, dashed line 
in Fig. These states are in the energy region of the 
Si PBS only, and show a small coupling with the Si con- 
tinuum. Inspection of the probability density of these 
states, shown in Fig. [2Jb) , indicates that they derive pre- 
dominantly from the atomic s-states of Si atoms in the 
second and first layer from the interface. The authors of 



□ Si □ GaAs 




FIG. 1: Dispersion of the interface states (thick lines) of the 
Si/GaAs (110) heterojunction relative to the PBS of the Si 
and GaAs bulk materials (hatched). The dashed line is a res- 
onant state localized on the first Si bilayer from the interface 
(see text). 



Ref. ^3 have found analogous resonant states in the case 
of the Ge/GaAs (110) interface. 

The band in the —7 to —4 eV energy range along the 
X-M line corresponds to states related to Si and Ga 
atoms at the interface (we called them Ssi-Ga states). 
The probability density of one of these states at M is 
displayed in two different planes, perpendicular and par- 
allel to the interface, in Figs.^e) and (e'), respectively. 
Such states contain a large s component, and also a non- 
negligible p contribution, of the Si and Ga atoms at the 
interface, and exhibit a strong antibonding character. 
They are pushed away from the Si-Ga bonding region 
into the more attractive Si-Si and Ga-As antibonding re- 
gions. We note that at the Ge/GaAs (110) interface, the 
corresponding states were found to have a dominant con- 
tribution from the Ga atom^ which is not the case here. 
We attribute this mainly to strain. Strain pushes the Si s 
states upwards in energy, giving rise to a stronger interac- 
tion with the Ga states, and hence to a more significant 
contribution of the group-IV atom to the antibonding- 
like interface state (this will become more clear from the 
study of the ultrathin Si (110) layers in the following 
sections) . 

The states of the next two bands that extend in the X - 
M-X' region (Psi-As and Psi-Ga) originate fromp states 
of the Si atoms at the interface and from p states of the 
interface As and Ga atoms, respectively (see Fig. |5Jc) 







FIG. 2: Contour plot of the probability density p (in 
e/supercell) in the (110) plane and planar average of the 
probability density (p) along the [110] direction of selected 
interface states of the Si/GaAs (110) heterojunction. (a) As- 
derived s state at X'. (b) Si-derived s state at X'. (c) Si and 
As derived p state at X '. (d) Si and Ga derived p state at X '. 
(e) Si and Ga derived s-p states at M . (e') The contour plot 
of the (e) state in the Si (110) plane adjacent to the interface 
(dashed line in the (e) plot). The probability density of the 
above states in atomic planes other than the ones considered 
is a small fraction of the one shown. 



and (d)), which form bonding orbitals across the inter- 
face. For the localization of the Psi-As states, the same 
considerations as for the Sa s states are valid. As far 
as the Psi-Ga states are concerned, we note that due to 
the repulsive ionic potential in the Si-Ga bonding region 
(with respect to the potential in the Ga-As or Si-Si re- 
gions), these states are pushed at the edge of (or at a 
higher energy than) the GaAs PBS. 
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FIG. 3: Contour plot of the probability density of two 
Si/GaAs (110) interface states in two different (110) planes 
containing Ga and As atoms (a) and Si atoms (b) at the in- 
terface, and corresponding planar average (p) along the [110] 
direction. Both states are at the X ' point of the zone. The 
dashed lines indicate the position of the planes where the con- 
tour plots were drawn. 



The remaining two types of state (Pca-As and Psi_si 
in Fig.^l are found only in the neighborhood of the X' 
point, and have their maximum probability density ei- 
ther on the GaAs (110) plane (Pca-As) or on the Si (110) 
plane (Psi-sO adjacent to the interface (see Fig.|3J). These 
states correspond essentially to As-Ga (Si-Si) p x bond- 
ing orbitals within the nearest-neighbor zigzag chains of 
the (110) plane. We note that the X' point corresponds 
to a zone-edge K vector perpendicular to the chains. In 
a tight-binding description^ this gives rise to a phase in 
the Bloch functions that leaves the sign of the p x orbitals 
of second nearest neighbor atoms in the chain unchanged, 
but reverses the sign of the atomic-like orbitals in adja- 
cent chains within the (110) plane. This leads to an anti- 
bonding configuration between the chains, but allows the 
formation of p x nearest-neighbor bonding orbitals within 
the chain, as seen in Fig. El Finally, the fact that the 
Psi-Si (Pca-As) states are lower (higher) in energy than 
the Psi-Ga (Psi-As) states can be explained by the more 
attractive ionic potential of Si as compared to Ga. 



IV. STATES INDUCED BY ONE SI (110) 
MONOLAYER IN GAAS 



Let us return to the Ssi states and also to the Ssi-Ga 
states of the Si/GaAs (110) heterojunction. Their for- 
mation and atomic character can be better understood 
if only few layers of Si are introduced in bulk GaAs. In- 
deed, the bulk PBS of Si is no longer present, which en- 
hances the localization of some of the Si induced states. 
We note that some of the states of the Si/GaAs inter- 
face are bound to be present also in the ultra-thin layer 
limit. In particular, all As and Ga-derived states can- 
not disappear since they essentially originate from GaAs 
band-edge features. We will see here and in the next sec- 
tion that actually most of the interface states of the fully 
developed Si/GaAs (110) heterojunction can be traced 
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□ GaAs 



1 Si ML 




FIG. 4: Dispersion of the localized states and resonances 
(thick lines) induced by 1 ML of Si (110) in GaAs relative to 
the GaAs bulk PBS. 



down to the ML limit. 

In Fig. 0] we have plotted the states induced by 1 ML 
of Si in the GaAs (110) homojunction, together with the 
GaAs PBS. The layer- induced states have been calcu- 
lated at k z — 0. However, the miniband width along 
k z is negligible given the thickness of the GaAs slab em- 
ployed in the present study. We have aligned the super- 
cell states with the bulk GaAs PBS using the calculated 
shift of the macroscopic average of the electrostatic po- 
tential in GaAs relative to its mean value in the super- 
cell. The group symmetry is C^, and the space group is 
symmorphic. Thus, all points and lines in the 2DBZ be- 
long to one-dimensional irreducible representations, and 
no additional degeneracy is introduced by time reversal 
symmetry. 

An Sas band is still present in Fig. 01 at approximately 
the same energy as in the interface case. This band de- 
rives from the s states of the As atoms adjacent to the Si 
layer that form symmetric states under the a xy reflection 
(see Fig.EIa)). The antisymmetric counterpart is found 
as a broad resonance (not shown in Fig. in the PBS of 
the bulk s valence band of GaAs. Conversely, a Si- related 
s band appears at about —9 eV at T , whose dispersion 
from r to X ' , to M , and to X is similar to that of 
the Ssi interface state. Again, this band corresponds to 
symmetric states with respect to <r xy . The probability 
density of one of these states (Ssi at X') is displayed 
in Fig. CJb). The probability density is predominantly 
on the Si atoms sitting on the anion sites in the con- 
tinuation of the zincblende lattice, but also reveals some 
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FIG. 5: Contour plot of the probability density p in 
the (110) plane and planar average of the probability den- 
sity (p) along the [110] direction of selected states of the 
GaAs/1 SiML/GaAs (110) heterostructure. (a) As-derived 
s state at X'. (b) Si-derived s state at X'. (c) Ga-derived 
s state at T. (d) Si and As derived p-state at X'. (e) Si 
and Ga derived p-state at X'. The probability density of the 
above states in atomic planes other than the ones considered 
is a small fraction of the one shown. 



bonding character with the nearest-neighbor Ga and Si 
atoms. Considering the Si layer as a perturbation to 
bulk GaAs^S these states separate from the GaAs s 
band and move up in energy, because of the repulsive 
on-site potential associated with the Si — > As substitu- 
tions. We note that because these localized states have 
their probability density mostly on the anion sites, within 
the Si (110) layer at the interface, they are qualitatively 
different from the Si bulk states present at similar en- 
ergies in Fig. ^ the latter states derive from Si atomic 
orbitals with identical weight on the two Si sublattice 
sites. The former states can thus give rise to Fano reso- 
nances when they interact with the continuum of bulk Si 
states, in the fully developed Si/GaAs (110) heterojunc- 
tion. The latter interaction generates bonding features 
between the Si atoms sitting on anion sites, in the layer 
closest to GaAs, and the Si atoms sitting on cation sites, 
in the neighboring layer on the Si side. This, in turn, 
pushes the probability density of the Ssi state towards 
the second Si layer, when the heterojunction is formed. 

In the neighborhood of the M point, along the T- 
X-M-X' path in Fig. |3| the Ssi state couples with the 
symmetric state we labeled Sc a - The latter state de- 
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FIG. 6: Left: contour plot of the probability density p 
in the (110) plane and planar average of the probability 
density (p) along the [110] direction of three states of the 
GaAs/1 Si ML/GaAs (110) heterostructure. Right: contour 
plot of the same states in the (110) plane indicated by the 
dashed line, (a) sp-Si and p-As derived state at X'. (b) Si 
and Ga-derived s state at M . (c) Ga-derived s state at M . 
The probability density of the above states in atomic planes 
other than the ones considered is a small fraction of the one 
shown. 



rive mainly from the s states of the Ga atoms adjacent 
to the Si layer (its probability density at T is shown in 
Fig. Et c ))- Consistent with the ordering expected from 
the corresponding s levels in the isolated atoms, the en- 
ergy positions of the localized s states of As, Si, and Ga 
are found in ascending order. The state labeled S' Ga , ob- 
served at M , is the antisymmetric counterpart of the Sc a 
state (see Fig. EI C )), while the symmetric Ssi-Ga state, 
found in the lower part of the GaAs stomach gap, has 
mixed Si and Ga s (with some p) atomic character (see 

Fig. EM). 

Another state (also symmetric under the a xy reflec- 
tion) is observed in the lower part of the GaAs stomach 
gap, but only near the X' point (SPgi-As)- The latter 
state is not present in the fully developed Si/GaAs (110) 
heterojunction, and originates mainly from p states of 
the As atoms adjacent to the Si layer and sp hybrids 
of the neighboring Si atoms (see Fig. Eta))- The Psi-As 
(Psi-Ga) state derives from the anti-symmetric combina- 
tion of the Psi-As (Psi-Ga) interface states on each side of 



the Si layer (see Figs. El and EJ. The Pca-As and Psi-Si 
states are symmetric states, which also derive from their 
interface counterparts (see Figs. and EJ. 

Finally, we find a well localized conduction state in 
part of the GaAs fundamental bandgap (Sg;_ As in Fig.QJ 
which does not exist in the fully developed Si/GaAs (110) 
heterojunction. The latter state is an antibonding Si- 
As state with strong sp 3 -hybrid character from the Si 
atoms occupying the cation sites within the (110) layer 
(see Fig. Etc)). This state is essentially a quantum-well 
state in its early stage formation. At larger Si intralayer 
thickness, the localization of this state, at specific K 
points of the 2DBZ, will be due to the conduction band 
offset of Si in GaAs (i.e., the quantum- well potential, 
within an effective- mass description). The Sg i As state is 
present within the GaAs bandgap in a large portion of 
the 2DBZ. Yet, we recall that the GaAs and Si bandgaps 
are affected by an error within the LDA-DFT approach, 
which may affect also the existence of the Sg iAs state. 
The calculated GaAs LDA (direct) gap is 0.9 eV, while 
the experimental valued is 1.52 eV (+0.62 eV correc- 
tion), and the Si LDA gap is about 0.5 eV to be com- 
pared to the experimental valued of 1.17 eV (+0.67 eV 
correction). The gap corrections, however, turn out to 
be similar for Si and GaAs; also, the corrections are ex- 
pected to be rather rigid shifts throughout the BZ. We 
therefore do not expect drastic changes in the alignment 
of the Si and GaAs conduction bands. To have an esti- 
mate of the conduction offset controlling the localization 
of the most bound Sg ; _ As states, we should compare the 
X (indirect) gap of GaAs (its calculated value is 1.4 eV; 
the experimental valued is 2.18 eV) with the correspond- 
ing gap of Si. Taking into account the 0.14 eV VBO, we 
thus conclude that the states near X' are bound by a 
conduction offset of about 0.87 eV, to be compared to 
the theoretical offset of 0.76 eV. This gives us confidence 
in the existence of the Sg;_ As states also in the actual 
quasiparticle interface spectrum. 



STATES INDUCED BY TWO SI (110) 
MONOLAYERS IN GAAS 



We end this study of non-polar (llO)-oriented struc- 
tures with the case of 2 MLs of Si in GaAs. This case is in- 
termediate between the fully developed interface and the 
single Si impurity-sheet perturbation. The point group 
symmetry of the 2 Si MLs in GaAs contains the reflection 
operations a yz , a xv , and a twofold rotation around the 
y axis. The point group is C% v , and the space group is 
non-symmorphic with fractional translation (i, ^,0) (in 
lattice units). Since our calculations do not take into 
account the non-scalar spin-orbit interaction, no addi- 
tional degeneracy is imposed along the T — X and V - X ' 
lines by time reversal symmetry. We thus expect along 
these lines pairs of symmetric and antisymmetric states 
for each atomic-orbital derived state. Conversely, the 
representations of the groups of X and X ' belong to the 




FIG. 7: Contour plot of the probability density p in the (110) plane indicated by the dashed line and planar average (p) along 
the [110] direction of three states of the GaAs/1 SiML/GaAs (110) heterostructure. All states are at the X' point of the zone. 
The dashed lines indicate the position of the planes where the contour plots have been drawn. 



two-dimensional irreducible representation X5, and that 
of M to the two-dimensional representation M 5 of C^v , 
for which an additional degeneracy is present. 

We plot in Fig. |H1 the states induced by 2 MLs of Si 
in bulk GaAs. We have aligned the localized states and 
resonances with the GaAs PBS only, as in the 1 ML case. 
In Fig. |H1 we distinguished the symmetric (by a xy reflec- 
tion) from the antisymmetric states, which have a small 
probability density between the adjacent Si layers. Com- 
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FIG. 8: Dispersion of the localized and resonant states in- 
duced by 2 MLs of Si (110) in GaAs relative to the GaAs 
bulk PBS. Solid lines are the symmetric (by a xy reflection) 
states. The dashed-dotted lines represent the corresponding 
antisymmetric states. 



paring the results for the 2 Si MLs with those obtained for 
1 Si ML and the fully developed interface, we see that the 
connection between the states of 1 ML and of the inter- 
face is made by the symmetric and antisymmetric states 
of the 2 MLs structure, which eventually become degen- 
erate for thicker Si layers. We note that the S'ca and 
SPsi-As states are no longer present already at 2 Si MLs 
coverage. Yet, we find a well localized Sg i As state in a rel- 
atively large portion of the 2DBZ. This state is localized 
in regions of the 2DBZ, which more closely correspond to 
the regions exhibiting a large Si/GaAs conduction band 
offset, as compared to the 1 ML case. Following the con- 
siderations given for the 1 ML case, this indicates that 
we arc in a situation in which the quantum well potential 
is nearly established, with a potential strength that is be- 
coming similar to the conduction-band offset between Si 
and GaAs. 



VI. STATES INDUCED BY A SI (100) BILAYER 
IN GAAS 

The situation is more complicated for (100)- than for 
(llO)-oriented heterovalent heterostructures, because of 
the polarity of the [100] growth direction. We mentioned 
in the Introduction that the fully developed abrupt 
Si/GaAs (100) heterojunction is thermodynamically un- 
stable, as it is macroscopically charged. Similarly, an 
abrupt Si (100) monolayer in GaAs is also macroscop- 
ically charged and thermodynamically unstable. An 
abrupt Si (100) bilayer, instead, gives rise to a dipolc, 
which is the lowest-order multipole which does not gen- 
erate a long-range electric field and is hence, in princi- 
ple, energetically acceptable. The dipole layer also corre- 
sponds to the highest-order multipole which can induce 
a change in the band alignment, and is therefore of possi- 
ble interest for interfacial band engineering. Here we thus 
focus on the bilayer case. Some caution, however, is in 
order when comparing to experiment. Although photoc- 
mission and internal photoemission studies have reported 
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some changes in the band alignment in the presence of 
Si interlayers; 7 i 8 i 9 this effect was found to saturate at Si 
coverages of about 0.5 ML, and to decrease in amplitude 
at higher coverages. This has been interpreted^ as an 
indication of significant deviations from a Si bilayer dis- 
tribution above 0.5 ML, due to atomic interdiffusion. It 
should be stressed, however, that the positions of the Si 
atoms are unknown in any case, and it is only specula- 
tively that one may assume a dipole layer distribution, 
even at low coverage. We thus consider the simplest pos- 
sible case of 2 MLs of Si in an abrupt bilayer configu- 
ration, and we will end this section with some consid- 
erations on the trends one may expect for bilayers with 
lower Si nominal concentration. Although 2 MLs of Si 
in a bilayer geometry is a structure which has not been 
produced experimentally so far, the study of this system 
is interesting in itself to establish a relationship between 
the interface states of the (110) and (100) heterostruc- 
tures. We also want to show that intrinsic scattering 
centers can exist in the (100) systems, even in an abrupt, 
defect-free heterostructure. 

In Fig. El we show the dispersion of the states induced 
by 2 MLs of Si in the (100) GaAs homojunction. The 
PBSs of the two GaAs bulks on each side of the junction 
(GaAs", GaAs d ) are rigidly shifted one with respect to 
the other by 1.15 eV, which is the calculated VBO of the 
relaxed heterostructure. In the LDA framework, the sys- 
tem is metallic near F . However, this would not be the 
case if the experimental gap of GaAs had been used. We 
note that the zone center has a very small weight in the 
BZ integrations (our grid actually does not include the 
r point). The charge transfer generated by the small 
valence-conduction band overlap at T has therefore a 
negligible influence on the band alignment (the computed 
potential lineup AV agrees within 0.1 eV with the pre- 
diction of the linear-response model of Ref.0). The point 
group symmetry of the isolated bilayer in GaAs is Ci v . 
Since we have two Si bilayers in the supercell, rotated by 
7r/2 one with respect to the other around the twofold z 
axis, a higher T)<id symmetry is present in the superlat- 
tice. The irreducible part of the 2DBZ of our calculation 
is thus one half of that of the isolated bilayer, and each 
point in the irreducible BZ includes states of both bi- 
layers. However no additional degeneracy is induced by 
time reversal symmetry. The bands are shown along the 
high symmetry lines connecting T = (0,0), J = {\,\) 
and K = (0,1) (in units of 2ir/a, where a is the GaAs 
lattice constant). 

Beginning from the lowest energies, we recover the As- 
related s-state band. In the present case, pairs of local- 
ized s states separated by AV appear, which correspond 
to As on the upper and lower (S As in Fig.lMa^ energy 
side of the heterostructure. We observe three different 
types of localized Si s-like states in the Si (100) bilayer 
structure: one derives mainly from Si atoms sitting on 
the Ga sites (Sg ; , see Fig. EJ C )), the second one is as- 
sociated mostly with Si atoms on the As sites (Sg ; , see 
Fig. Eld)), while the third one involves Si atoms from 
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FIG. 9: Dispersion of the localized and resonant states (thick 
lines) induced by 2 MLs of Si in the (100) GaAs homojunc- 
tion relative to the PBS of the GaAs bulks, GaAs" and 
GaAs d , which are shifted one with respect to the other by 
the dipole potential induced by the Si interlayer. The dashed 
line indicates the state corresponding to the Sg» state of the 
GaAs/1 SiML/GaAs (110) homojunction (see text). 



both cation and anion sites (Ssi-Si m Fig. HOf bj j. and 
shows strong Si-Si bonding features. The first two types 
of states are found essentially around the zone edge K 
point, whereas the Ssi-si states are observed in the re- 
maining part of the 2DBZ. The character of the Sgi-Si 
states changes progressively along the J-K direction, 
and matches the character of the Sg ; and Sgj states at 
K. 

The Ga-related localized s state of the (110) 
GaAs/1 Si ML/GaAs heterostructure is found as a res- 
onance in the present case (dashed line in Fig.^J). Con- 
versely, we find well-localized states which derive from 
interacting Ga s and Si p states (SPca-Si in Fig. IllT ajj. 
We also find an SPAs-Si state, which results from the in- 
teraction between As sp hybrids and Si p orbitals (SPAs-Si 
in Fig. rnT bil. The in-plane (110) Pca-As interface state 
is no longer present, as one might have expected since 
the planes of Ga-As zigzag chains are now perpendicu- 
lar to the interface. The Psi-Si state, instead, is found 
as a bilayer state, oriented in the (100) direction (see 
Fig. I12f bj j. We note that the Psi-Si state appears only 
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FIG. 10: Contour plot of the probability density p in the 
(110) plane and planar average (p) along the [100] direction 
of some states in the GaAs/2 Si MLs/GaAs (100) heterostruc- 
ture. (a) GaAs d As-derived s-state at J . (b) Si-Si-derived 
s-state with lowest energy at J . (c) Si d -derived s-state at K . 
(d) Si u -derived s-state at K . The probability density for the 
above states in atomic planes other than the ones considered 
is a small fraction of the one shown. 



near the valence edge of the GaAs d PBS, along the T- J 
line. The state originating from the interaction between 
the p states of As and Si (PAs-Si in Fig. IT2T a)) is found 
near the zone-edge K point, and is a resonance in the 
GaAs d PBS. 

Finally, we also report two states which are present in 
the main PBS gap (apart from the metallic region near 
r ) . The first one is a Ga-Si bonding state resulting from 
the interaction between the p states of Ga and Si (PGa-Si> 
see Fig. ll2f c)) at the GaAs"/Si junction. The second one 
is a Si- As antibonding state deriving from sp 3 hybrids of 
the Si and (to a lesser extent) of the As atoms (S Si _ As , see 
Fig. I12f d)) at the Si/GaAs d junction. Both states, and 
in particular the Pea-Si state, exhibit a relatively weak 
dispersion along the T-J line, leading to an enhanced 
DOS near the GaAs u valence-band edge (P Ga-Si) and the 
GaAs d conduction-band edge (Sg;_ As ). 

In the presence of only one of the two GaAs u /Si and 
Si/GaAs d interfaces, we expect the resulting Pea-Si or 
Sgj _ As states to be within the GaAs optical gap near 
r . This is confirmed by ab initio computations we per- 
formed for related engineered systems, namely abrupt 
Al/6 Si MLs/GaAs (100) metal-semiconductor structures 
with either an As- or a Ga-GaAs (100) termination at 
the polar Si/GaAs interfaced The dispersion of the cor- 
responding S Si _ As and Pea-Si states is shown in Fig. ^2 
along the high symmetry direction of the 2DBZ, together 
with the probability density of these states. Although 
the As- and Ga-terminated Al/6 Si MLs/GaAs (100) 
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FIG. 11: Left: contour plot of the probability density p in the 
(110) plane and planar average (p) along the [100] direction of 
two states in the (100) GaAs/2 Si MLs/GaAs heterostructure. 
Right: contour plot of the same states in the (100) plane 
indicated by the dashed line, (a) s-Ga and p-Si derived state 
at J . (b) sp As and Si-derived state at K . The probability 
density for the above states in atomic planes other than the 
ones considered is a small fraction of the one shown. 



metal/semiconductor structures incorporate only one of 
the two types of charged Si/GaAs (100) interfaces each, 
they exhibit no macroscopic electric field in the bulk 
semiconductor (GaAs) — because of the presence of the 
metal — and display either a GaAs d energy-shifted bulk 
(in the As-terminated junction) or a GaAs" energy- 
shifted bulk (in the Ga-terminated junction). Also, the 
thickness of the Si (6 Si monolayers pseudomorphically 
strained to GaAs) insures negligible interaction between 
the continuum of Al bulk states, on one side of the inter- 
layer, and the resonant Sg;_ As or Pea-Si Si/GaAs-interface 
states, on the other side of the interlayer. In these condi- 
tions, the As-terminated (Ga-terminated) junction gives 
rise (see Fig.[I3J> to a band of the S Si _ As (Pea-Si) Si/GaAs 
interface states which has a dispersion similar to that of 
Fig. within the GaAs PBSs bandgap, but with an en- 
ergy at r which is well within the gap, namely ~ 0.2 eV 
below (~ 0.5 eV above) the conduction- (valence-) band 
edge of the GaAs PBS. 

The Sg;_ As and Pea-Si states that we find at the abrupt 
As-terminated and Ga-terminated Si/GaAs (100) inter- 
faces, respectively, are the (100) analogues of the Sg;_ As 
and Psi-Ga states that we observed in the non-polar 
GaAs/Si/GaAs (110) heterostructures. However, con- 
trary to the case of the globally neutral Si/GaAs (110) 
junction, at the As- terminated Si/GaAs (100) inter- 
face, the attractive potential created by the globally 
positive excess ionic charge of the As-Si donor bonds 
is sufficiently strong to pull down Sg;_ As conduction 
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FIG. 12: Contour plot of the probability density p in the 
(110) plane and planar average (p) along the [100] direction 
of some states in the (100) GaAs/2 Si MLs/GaAs heterostruc- 
ture. (a) GaAs d As and Si d -derived p-state at J . (b) Si- 
derived p-state at J . (c) GaAs" Ga and Si-derived p-state at 
J . (d) Si d -derived antibonding sp-state at J . The probability 
density for the above states in atomic planes other than the 
ones considered is a small fraction of the one shown. 



states within the GaAs fundamental gap. Similarly, 
the potential created by the globally negative ionic 
charge of the Ga-Si acceptor bonds is able to push 
U P Pea-Si states from the valence band into the GaAs 
gap. Hence, the Sg ; _ As (Pea-Si) states become poten- 
tial scattering centers at the polar As-terminated (Ga- 
terminated) Si/GaAs (100) heterojunction. 

As mentioned before, Si may be confined in a bilayer 
geometry only at nominal coverages lower than or equal 
to 0.5 ML. We have therefore also investigated what hap- 
pens to the Sg;_ As and PGa-Si states when the Si coverage 
is reduced to 0.5 ML. We have considered two different 
ordered GaAs/Si/GaAs (2 x 2) interlayer atomic config- 
urations corresponding each to one Si — > Ga and one Si 
— > As atomic substitution per (2 x 2) surface-unit cell 
within the GaAs bilayer. In the first configuration, the 
pair of Si atoms sit on adjacent Ga-As sites, so that each 
Si atom has three As or three Ga nearest neighbors and 
one Si nearest neighbor (configuration A), whereas in the 
second configuration each Si atom has four As or four Ga 
nearest neighbors, and no Si nearest neighbor (configu- 
ration B)^4 

In both A and B configurations, wc find that the Sg i As 
(PGa-Si) states associated with the Si- As (Si-Ga) donor 
(acceptor) bonds become resonances at L. The energy 
positions of the resonances depend somewhat on the Si 
arrangement. The Sg i As resonance is located at +0.2 eV 
from the conduction band edge in the A configuration, 
and at +0.1 eV from the same edge in the B configu- 



ration. The PGa-Si resonance remains, instead, a near- 
band edge resonance (within 0.1 eV from the valence 
band edge) in both cases. We note that, in both cases, 
the conduction resonance gives rise to a stronger DOS 
feature than the valence resonance, due to the lower 
density of bulk states present near the conduction edge. 
The conduction resonance may be responsible for a pe- 
culiar feature observed in internal photoemission experi- 
ments by dell'Orto et al~ The associated increased DOS 
at L is indeed detectable, in principle, by transport 
measurements ^ 



VII. KOSTER-SLATER MODEL FOR THE 
INTERFACE-STATE PROBLEM 

The ab initio calculations of the preceding sections 
have provided us with the energy position and atomic- 
like character of various states localized at the Si/GaAs 
interface. However, it is not simple to extract from the 
full-fledged ab initio calculations a general understand- 
ing of the mechanisms of interface-state formation. For 
example, it is evident from Fig. ^ that except for the 
S As states which form a complete band along the high- 
symmetry lines, the other states appear only along spe- 
cific lines or near specific points of the 2DBZ. This is 
a common feature observed for many different interfaces 
studied in the literature. In this section, we thus try to 
rationalize the above results based on a more general, yet 
simplified description, that contains the main physical 
ingredients of the problem. Specifically, one would like 
such a model to include the effects of (?) the strength and 
sign of the local potential at the interface generated by 
the different chemical nature of the interface atomic con- 
stituents, (ii) the band structures of the two bulk materi- 
als, and (Hi) the macroscopic lineup of the two bulk band 
structures across the interface. We report below a simple 
analytical criterion for the existence of interface states, 
which explicitly relates the existence of these states to 
(i), (ii), and (Hi). A derivation of the analytical form of 
this criterion, which is mathematically simple and phys- 
ically transparent, is also presented. A detailed and rig- 
orous quantum-mechanical derivation, together with the 
explicit expressions for the microscopic interface-specific 
terms, is given in Ref. fist 

An interface can be regarded as a perturbation of 
the bonding structures of two different bulk materials. 13 
Also, since for stable interfaces the perturbation involved 
is neutral and short ranged, we may formulate the prob- 
lem in terms of few perturbed layers, and thus derive a 
criterion similar to the one Koster and Slater introduced 
for deep impurities in semiconductors i^i Following the 
authors of Ref. fl3l we construct the interface as a per- 
turbation of two bulk semiconductors A and B. We write 
the unperturbed Hamiltonian as 
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FIG. 13: Resonant Si/GaAs-interface states in the abrupt Al/6 SiMLs/GaAs (100) junctions with As- (left) and Ga- (right) 
terminated GaAs (100) at the Si/GaAs (100) interface, (a) Dispersion of the interface states along the high symmetry directions 
of the 2DBZ (solid line); the gray area corresponds to the projected bulk band structure of GaAs. The zero of energy corresponds 
to the GaAs valence band maximum, (b) Probability density of the resonant states at the F point. The contour plots are 
shown in (010) planes including the As-Si-As (left) and Ga-Si-Ga (right) bonds at the Si/GaAs interface. Solid (dotted) 
circles and lines indicate atoms and bonds which are in (out) of the plane. The planar average of the density is shown below 
the contour plots. 



where Hq and Hq are the Hamiltonians of the two bulk 
periodic crystals, with standard periodic boundary con- 
ditions over N layers. The corresponding one-particle 
Green's function of the system is 



G = 



G£ 
Gq 



(2) 



with Gq and Gq the Green's functions of the Hamiltoni- 
ans Hq and Hq , respectively. A perturbation U (which 
will be specified below) is then introduced; U acts on this 
system, and in particular couples A and B, and yields an 
interface A/B (plus two free surfaces). The total Hamil- 
tonian describing the perturbed system is H = Hq + U. 
A basis set of Wannier functions is assumed in the follow- 
ing discussion. Specifically, as U preserves the periodicity 
in planes parallel to the interface, the basis functions at 
a given K -point in the 2DBZ consists of localized-layer 
Wannier functions^ in A and B, obtained from the Bloch 
functions of the corresponding bulk with wave- vector pro- 
jection in the 2DBZ equal to K . 

The construction of the interface, and thus of the per- 
turbation U, is schematized in Fig. El We consider the 
Si/GaAs (110) interface, and assume for simplicity that 
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FIG. 14: Schematic description of the interface creation pro- 
cedure discussed in the text. 



the interactions are described by tight-binding nearest- 
neighbor Hamiltonians. To keep the expressions as sim- 
ple as possible, we also refer in the following to one band 
£ A(B){K,k z ) for each bulk semiconductor. The gener- 
alization to the multiband case is straightforward^ A 
layer of atoms of material A (e.g., GaAs) and a layer of 
atoms B (Si) are first removed on each side of the inter- 
face in bulk space A and B, respectively. This procedure 



12 



is analogous to the creation of free surfaces and can be 
performed in practice by changing the on-site parame- 
ters Va and Vb , which eventually are let go to infinity*^ 
Then, the interactions between the atoms of the left-hand 
side A and right-hand side B surfaces (Si-Ga and Si- 
As) are switched on across the interface. These interac- 
tions are represented by off-diagonal terms Uab, which 
are nearest-neighbor tight-binding matrix elements cou- 
pling the two interface layers. Finally, the charge trans- 
fer at the interface, i.e., the perturbation induced by the 
interface relative to a superposition of bulk crystal po- 
tentials, is added as on-site diagonal matrix elements Ua 
and Ub for the semiconductor interface layer A and B, re- 
spectively. Since these perturbations are measured with 
respect to the average electrostatic potential of the cor- 
responding bulks, the lineup is introduced in the Green's 
functions of the bulks. If we assume that the removal 
procedure described above has been performed, the only 
relevant block of the matrix U is 

U=( U , A UAB ), (3) 

where all elements are scalars in the one-band model. 
The condition for having bound states can then be de- 
rived from the determinantal equation^ 

det[l - G (E)U] =0, (4) 

and reads 

1 - U A G A - U B G% - (\U AB \ 2 - U A U B )G A G® = 0. (5) 

If the zero of energy is fixed at the top of the valence 
band EA(K,k z ), the Green's functions are: 

G o = Jf E E ^ kz y (6) 

and 

C B = — X~^ I (7) 

N ^ E-e B (K,k z )-VBO(A/B)- {) 

The condition (JSJ) tends to the usual Koster-Slater re- 
lation, 1 — UGo = 0, for deep impurities (U being the 
potential perturbation induced by the impurity in the 
host described by Go) when only one semiconductor and 
on-site interactions are retained. Given the linear rela- 
tion, Gq(E) — J deD(e)/(E — e), between the Green's 
function of a bulk and the corresponding electronic den- 
sity of states (DOS), D(e), we find from Eq. (|SJ) that, for 
a given set of parameters (Ua, Ub and Uab), the exis- 
tence of bound states and resonances depends on: (?) the 
strength of the DOS of each bulk band with respect to the 
corresponding on-site parameter U^ B y and (ii) the am- 
plitude of the product of the DOS's of the two bulk bands 
with respect to the interface coupling parameter |C/4b|~ 2 - 
Localized states may appear above/below the bulk bands 



according not only to the strength of the on-site terms Ua 
and Ub at the interface, but also to the relative strength 
of these terms and the coupling terms Uab*^ In fact, 
when Ua = Ub = 0, Uab can push states both above 
and below the bulk bands. It should be noted, however, 
that when the band offset is non-vanishing, even in a 
one-dimensional model, the on-site perturbations gener- 
ated by the interface may induce resonances rather than 
bound states. This is in contrast to the case of deep im- 
purities, where a state is always bound in one dimension 
(for local potentials). 

With the help of Eq. ©, we can now try to identify 
the regions and points of the 2DBZ where interface states 
are most likely to appear. In view of the relationship 
between Green's function and DOS, we have examined 
the DOS of both Si and GaAs bulks, for a given point 
of the (110) 2DBZ, integrated along a line in fc-spacc 
parallel to the growth direction. The result is plotted in 
Fig- 113 for the four high-symmetry points of the (110) 
2DBZ. The DOS has been calculated using tetragonal 
(properly strained) supercells of 4 atoms oriented along 
the (110) direction for both bulks. For each point in the 
2DBZ we sampled the (110) direction with 26 fc-points 
distributed on a uniform grid. A Gaussian broadening of 
0.1 eV was used. We have employed the same procedure 
as in the previous section to align the DOS's of the two 
bulks. 

The s DOS of GaAs is a sharp feature, well separated 
from the Si bands at X, M, and X'. The localized Sa s 
state is induced thus by an on-site perturbation. The 
Uas on-site potential is attractive? 1 ^ due to the attrac- 
tive potential of the nearest Si atoms, and is responsible 
for the existence of this state. The Ssi_Ga states origi- 
nate mainly from both the Si and GaAs peaks located at 
about —6 eV and —7 eV in the X and M DOS's spectra, 
respectively. These edges are associated with Si (GaAs) 
electronic states that contain an important s component, 
together with a non-negligible p contribution, from the 
Si (Ga) atoms. The on-site Uq& and Usi matrix elements 
of adjacent Ga and Si atoms are repulsive^ However, 
contrary to the case of the As s states, the Si and GaAs 
edges are separated by less than 1 eV, and the Uab term 
also contributes to the formation of the Ssi-Ga states. At 
r and X', instead, no comparable DOS features are de- 
tected at similar energies. 

At X', two pronounced and strongly overlapping Si 
and GaAs DOS features are visible: with band edges 
at about —2.5 eV and —1 eV. The lower edges con- 
tribute predominantly to the Psi-As and Pca-As states, 
and the upper ones to the Psi-Si and Psi-Ga states. How- 
ever, only the Psi-As and Psi-Ga states can be followed 
along the 2DBZ up to the X point. The strong overlap 
near the band edges make the contribution of the Uab 
term important, as evidenced by the mixed Si-GaAs na- 
ture of the resulting interface states. These states are 
localized by the attractive (repulsive) ionic potential in 
the Si-As (Si-Ga) bonding region^ At the T point, no 
bound state is present in the fundamental gap. This is 
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FIG. 15: Wave-vector resolved density of states (DOS) of bulk Si and GaAs for the four high-symmetry points of the (110) 
2DBZ, and integrated along a line in fc-space parallel to the [110] crystallographic direction. The two DOS are aligned using 
the calculated VBO of 0.14 eV at the Si/GaAs (110) interface. The volume of normalization is the GaAs unit-cell volume. The 
degeneracies of the Si (GaAs) valence-band DOS features, given in order of increasing energy, are: 1 (1, 1) at F , 2, 2 (1, 1, 2) 
at X and M, and 1, 1, 2 (1, 1, 2) at X'. 



due to the attractive character of the Ua s and Usi ma- 
trix elements^ and to the finite value of the band off- 
set. The off-diagonal term U ab has thus not enough 
strength to push a state in the gap; we note that the sit- 
uation changes, instead, at the polar, negatively charged 
Ga-terminated Si/GaAs (100) interface, where the re- 
pulsive potential induced by the negative charge pushes 
the Psi-Ga state within the GaAs bandgap. Similar con- 
siderations apply to other abrupt interfaces studied in 
the Iiterature?ii*i2iiiii2ii£ii§i2! giving us confidence in the 
soundness of the above model description. We also note 
that estimates for the interface-bonding parameters may 
be obtained from relatively simple models, such as the 
tight-binding approach of Ref. [l3l 



VIII. SUMMARY AND CONCLUSIONS 

Using an ab initio pseudopotential approach, we have 
studied the states induced by thin Si interlayers in 
abrupt GaAs/Si/GaAs (110) and (100) heterostructures. 
We have also investigated the interface states of the 
fully developed Si/GaAs (110) heterojunction and reso- 
nant Si/GaAs (100) interface states occurring within the 
GaAs bandgap in Al/6 SiMLs/GaAs (100) heterojunc- 
tions. We have examined the bonding properties and 



atomic character of these various states, and their possi- 
ble connections. 

The results reveal interesting common features be- 
tween localized states and resonances in structures with 
different interface orientations and, for a given orienta- 
tion, between states localized by 1 or 2 Si monolayers and 
the interface states of the fully developed Si/GaAs junc- 
tion. By studying the evolution of the resonances and lo- 
calized states with Si-interlayer thickness, we were able to 
explain the origin of Si-related resonant interface states 
occurring in the fully developed Si/GaAs (110) junction. 
The ab initio results also indicated that most of the lo- 
calized states of the Si/GaAs (110) interface persist up 
to 2 and even 1 monolayer Si-interlayer coverage. More- 
over, the most striking valence- (conduction-) near-band 
edge states found in the GaAs/2 SiMLs/GaAs (100) 
heterostructure were also shown to originate from lo- 
calized As-Si donor-bond (Ga-Si acceptor-bond) states 
of the fully developed As-terminated (Ga-terminated) 
Si/GaAs (100) interface. 

In order to gain further insight into the mechanism of 
formation of interface states, we have presented a Koster- 
Slater type model developed for the interface-state prob- 
lem. A simple analytical criterion for the existence of 
interface states [Eq. JSJ] has been derived, which allows 
us to predict the existence and/or energy of interface 
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states from the DOS structure of the constituent materi- 
als and from some interface bonding parameters^ With 
this model, we were able to explain the different type of 
localized states obtained from the ab initio calculations 
for the GaAs/Si interface, and account for their energy 
position and location within the Brillouin zone. More 
generally, this model is expected to be useful to predict 
trends due to changes in the interface atomic structure 
(e.g., due to atomic intermixing^) or in the electronic 
band structure of the bulk materials (e.g., due to alloy- 



ing or different polytypes^I). 



Acknowledgments 

We acknowledge the hospitality and support of the 
Ecole Polytechnique Federale of Lausanne where most 
of this work has been performed while two of us (MD 
and CB) were completing their PhD thesis. 



1 See e.g., A. Franciosi and C. G. Van de Walle, Surf. Sci. 
Rep. 25, 1 (1996); and references therein. 

2 F. Capasso, Science 235, 172 (1987). 

3 W. A. Harrison, E. A. Kraut, J. R. Waldrop, and R. W. 
Grant, Phys. Rev. B 18, 4402 (1978). 

4 M. Peressi, S. Baroni, R. Resta, and A. Baldereschi, Phys. 
Rev. B 43, 7347 (1991). 

5 S. De Franceschi, F. Altomare, F. Beltram, M. Lazzarino, 
L. Sorba, and A. Franciosi, J. Appl. Phys. 84, 4637 (1998). 

6 M. Moreno, H. Yang, M. Horicke, M. Alonso, J. A. Martm- 
Gago, R. Hey, K. Horn, J. L. Sacedon, K. H. Ploog, Phys. 
Rev. B 57, 12314 (1998); M. Moreno, J. L. Sacedon, M. 
Alonso, M. Horicke, R. Hey, J. Avila, M. C. Asensio, K. 
Horn, K. H. Ploog, Phys. Rev. B 58, 13767 (1998). 

7 L. Sorba, G. Bratina, A. Antonini, J. F. Walker, M. 
Mikovic, G. Ceccone, and A. Franciosi, Phys. Rev. B 43, 
2450 (1991); G. Bratina, L. Sorba, A. Antonini, G. Biasiol, 
and A. Franciosi, Phys. Rev. B 45, 4528 (1992). 

8 T. delPOrto, J. Almeida, C. Coluzza, A. Baldereschi, G. 
Margaritondo, M. Cantile, S. Yildirim, L. Sorba, and A. 
Franciosi, Appl. Phys. Lett. 64, 2111 (1994). 

9 M. Marsi, R. Houdre, A. Rudra, M. Ilegems, F. Gozzo, 
C. Coluzza, and G. Margaritondo, Phys. Rev. B 47, 6455 
(1993). 

10 S. Modesti, D. Furlanetto, M. Piccin. S. Rubini, and A. 
Franciosi, Appl. Phys. Lett. 82, 1932 (2003). 

11 G. A. Baraff, J. A. Appelbaum, and D. R. Hamann, Phys. 
Rev. Lett. 38, 237 (1977); G. A. Baraff, J. A. Appelbaum, 
and D. R. Hamann, J. Vac. Sci. Technol. 14, 999 (1977). 

12 W. E. Pickett and M. L. Cohen, Phys. Rev. B 18, 939 
(1978); W. E. Pickett, S. G. Louie, and M. L. Cohen, Phys. 
Rev. B 17, 815 (1978). 

13 J. Pollmann and S. T. Pantelides, Phys. Rev. B 21 709 
(1980). 

14 T. Saito and T. Ikoma, Phys. Rev. B 45, 1762 (1992). 

15 L. Zhu, E. Wang, and L. Zhang, Phys. Rev. B 56, 10308 
(1997); and references therein. 

16 A. Laref, H. Aourag, B. Belgoumene, and A. Tadjer, J. 
Appl. Phys. 94, 5027 (2003). 

17 M. Peressi, F. Favot, G. Cangiani, and A. Baldereschi, 
Appl. Phys. Lett. 81, 5171 (2002); and references therein. 

18 M. Di Ventra, C. Berthod, and N. Binggeli, Phys. Rev. B 
62, R10622 (2000). 

19 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 



(1980). 

20 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

21 L. Kleinman and D. M. Bylander, Phys Rev Lett. 48, 1425 
(1982). 

22 H. J. Monkhorst and J. P. Pack, Phys. Rev. B 13, 5188 
(1976). 

23 Numerical Data and Functional Relationships in Science 
and Technology, Landolt-Bornstein, New Series, Vol. 22 
pt. a (Springer, 1987). 

24 See, e.g., M. Peressi and S. Baroni, Phys. Rev. B 49, 7490 
(1994). 

25 C. G. Van de Walle, Phys. Rev. B 39, 1871 (1989). 

26 We note that localized interface states may exist also 
within the bulk PBS if their symmetry differs from that 
of the bulk states. This is, however, relatively unlikely 
in semiconductor heterojunctions, especially at generic k 
points of the BZ. 

27 A. Baldereschi, M. Peressi, S. Baroni, and R. Resta, in 
Proceedings of the International School of Physics "Enrico 
Fermi" Course CXVII, Semiconductor Super-lattices and 
Interfaces, edited by L. Miglio and A. Stella (Academic, 
New York, 1993), p. 59. 

28 M. Peressi, N. Binggeli, and A. Baldereschi, J. Phys. D 31, 
1273 (1998). 

29 G. F. Koster, J. O. Dimmock, R. G. Wheeler, and H. Statz, 
Properties of the Thirty-two Point Groups (MIT Press, 
Cambridge, 1963). 

30 W. A. Harrison, Electronic Structure and the Properties 
of Solids (W. H. Freeman and Company, San Francisco, 
1980), ch. 3. 

31 G. F. Koster and J. C. Slater , Phys. Rev. 96, 1208 (1954). 

32 M. Di Ventra and K. A. Mader, Phys. Rev. B 55, 13148 
(1997). 

33 C. Berthod, Ph.D. thesis, EPFL (1999). 

34 Configuration A is found to be lower in energy by ~ 0.2 eV 
per Si atom than configuration B. 

35 T. dell'Orto, M. Di Ventra, J. Almeida, C. Coluzza, and 
G. Margaritondo, Phys. Rev. B, 52 R2265, (1995). 

36 E. N. Economou, Green's Functions in Quantum Physics, 
(Springer, 1979). 

37 M. Di Ventra, Appl. Phys. Lett. 79, 2402 (2001). 



